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Abstract 

High resolution simulations reveal that in the cold dark matter scenario the structures form 
hierarchically and a large number of substructures survive in the galactic halos. The substructures 
can be probed if they emit gamma rays via dark matter annihilation. We calculated the gamma ray 
fluxes from the dark matter annihilations in the substructures of our Galaxy within the frame of the 
minimal supersymmetric extension of the standard model. The uncertainties of the prediction from 
both the low energy supersymmetry and especially from the density profiles of dark matter in the 
substructures are carefully investigated. The cumulative number of substructures emitting gamma 
rays above any given flux is calculated. Detectability of the gamma rays from the substructures is 
discussed. We propose the viability to detect these signals through the ground large field of view 
detectors. 



I. INTRODUCTION 

Present observations strongly support a standard cosmology (ACDM) in which the 
Universe is spatially fiat and its energy budget is balanced by ~ 4% baryonic matter, ~ 23% 
cold, coUisionless non-baryonic dark matter (CDM), and ~ 73% dark energy. In the ACDM 
cosmology, the luminous galaxies and clusters of galaxies form within the halos of CDM 
where the dark matter potential wells trap the baryonic gas, which eventually cools and 
condense to form the galaxies. The DM halos are assumed to form hierarchically bottom 
up via gravitational amplification of initial fluctuations, generated during an early epoch of 
inflation with a nearly scale invariant primordial power spectrum. In this paradigm, small 
mass objects collapse first and merge into larger and larger halos over time. 

When a small halo merges into a large host system it is not immediately destroyed, but 
instead begins to orbit within the host, gradually losing mass to the parent halo due to 
the action of the tidal force from the host, the dynamical friction and the close encounters 
with other subhalos. N-body simulation has been extensively used to study the merging 
history and the structure of the CDM halos. As the recent development of fast algorithms 
to integrate the orbits of millions of particles, high resolution simulations indicate that a 
fraction of about 10% of the total halo mass may have survived tidal disruption and appear 
as distinct and self-bound substructures or subhalos inside the virialized host halos 3, \4, nl 
BBBBQ. The «.a. confirmation. a.e self ..ni.a. with a smooth host halo and a s.aU 
fraction of the total mass in subhalos. The halos that contain a wealth of substructures 
resemble the observed clusters which host the galaxies, while the internal structure of a 
galaxy-sized halo that hosts galaxy satellites looks just like a rescaled version of a rich 
cluster. 

The existence of substructures in the halos have been confirmed by different high reso- 
lution numerical simulations as a generic picture of the CDM cosmology with hierarchical 
structure formation. However, despite the great success of the CDM scenario in describing 
both the large scale distribution of matter in our Universe and the structure of galaxies 



and clusters 



12| |. simulation shows that halos similar to that of Milky Way (MW) should 



host hundreds of subhalos, which apparently overpredicts the abundance of substructures 
by an order of magnitude compared with the 11 observed dwarf galactic satellites of the 
MW This CDM problem on sub-galactic scale is regarded as one of the most 
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fundamental issues that has to be addressed. 

Several possible resolutions have been proposed to this apparent discrepancy. One pro- 

: c matter 
|l7| . and 

non-thermally produced dark matter modeljisl. Q|. Another possibility is to feature the 



posal is to change the nature of the dark matter, including self-interacting dar 
modelQ], warm dark matter model IigI. annihilating dark matter model 



inflation potentia 
subgalactic halos 



that suppress small scale power and thus reduce the predicted number of 



3, 



21[. However, these models become gradually disfavored by the recent 



observations and numerical simulations 22, |23[ . The astrophysical mechanisms explain the 
discrepancy by suppressing dwarf galaxies formation in subgalactic halos and claim that 
only very massive substructures contain stars and most substructures are dark. Therefore 
detection of the non-luminous subhalos in the galactic halos through lensing effects j25|. 
tidal streams 2^, or any other method promise to distinguish between these alternatives. 
It seems that to account for the flux anomalies observed in radio lensing the amount of 
substructures predicted by the ACDM model is required {2^. 

The subhalos may also be lit up by the annihilation of DM into 7-rays and probed by 7-ray 
detectors if the DM particles are in form of weakly interacting particles [28]. The subhalos 
in MW greatly enhance the fluxes of the annihilation products and therefore enables us to 
detect these products, since they produce many denser regions in the smooth background. 
To predict the intensity of the 7-ray fluxes it is necessary to study the nature of the DM 
particles. 

From the point of view of particle physics, the existence of non-baryonic dark mat- 
ter clearly indicates the new physics beyond the standard model (SM) of particle physics. 
Among the large amount of candidates proposed for non-baryonic DM the leading scenario 
involves the weakly interacting massive particles (WIMPs), which is well motivated by the- 
oretical extensions of the SM. The weakly interacting relics from the early Universe with 
the WIMP mass from some tens of GeV to several TeV can naturally give rise to relic den- 
sities in the range of the observed DM density. The most popular and natural extension of 
the SM seems to be its supersymmetric (SUSY) version, i.e., the minimal supersymmetric 
standard model (MSSM). The lightest supersymmetric particle (LSP) of the MSSM, usually 
the neutralino, is neutral and stable due to R-parity conservation and provides an excellent 
candidate of CDM. Search for dark matter via detection of its annihilation secondaries is 
strongly motivated to unveil the form of new physics beyond the SM. 
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Assuming that neutralino makes up the dark matter, in this paper we calculate the 7-ray 
fluxes from the dark matter annihilations in the MW subhalos. There have been similar 
studies in the literature 
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3ai3ii. 

In this work we have addressed the uncertainties from 



the unknown low energy SUSY parameters and from the not well determined DM density 
profile in subhalos simultaneously. Especially we pay more attention on how to determine the 
distribution and the density profile of the subhalos. Our results show that the prediction 
depends heavily on both sides of particle physics and the cosmological evolution of the 
large scale structure. The constraints on the DM annihilations from the observations of 
EGREtQ, CANGAROO 110 and HESS ^\ are taken into account. The possibility of 
detection of these 7-rays is then discussed. Complementary to the space experiments, such 



as GLAST |35|| and the atmosphere Cerenkov detectors, such as VERITAS l36ll . MAGIC l3'i 
or IIESS|38j, we find the large ground cosmic ray arrays, such as the ARGOj39| and the 
IIAWCj40| project may have the ability to detect the annihilations from the subhalos if the 
neutralino is as heavy as about 500 GeV and the density profile of subhalos has a steep 
central cusp. Therefore in most calculations we specify the quantities corresponding to the 
ground array detectors, such as the angular resolution and the threshold energy. 

The paper is organized as follows. In the Section II, we present the general formulas for 
calculating the 7-ray fluxes from DM annihilations. In the Section III, we calculate the 7-ray 
flux from the DM annihilations at the Galactic Center (GC) and study the uncertainties 
from the 'particle factor' of Eq. (j2)). Since there are several observations at the GC, we 
consider the constraints on the SUSY parameters from these observations. In Section IV 
we present our results of the 7-ray fluxes from subhalos after determining the distribution 
and density profile of the subhalos. The detectability of the signals by different types of 
detectors is discussed in Section V. In Section VI we give our conclusions. 



II. FLUXES BY DARK MATTER ANNIHILATION 

It is easy to show that the average annihilation rate in unit time and unit volume is given 

by 

flH->nV2^<^ (1) 

where a and v are the annihilation cross section and the relative velocity of the two dark 
matter particles respectively, n and p are the number and mass densities of dark matter and 
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m is its mass. We note that the annihilation rate is proportional to the square of the dark 
matter density and therefore, a high density region can greatly enhances the annihilation 
fluxes. 

The radiation fluxes from a dark matter halo is therefore given by 

in unit of 1 particle GeV^^cm^'^s^^, where d is the distance from the detector to the source 
where dark matter annihilates and (f){E) is the differential flux at energy ii^ by a single 
annihilation in unit of 1 particle GeV~^. At the last step of Eq. the integration 

is given alone the line-of-sight I, which is related with the galactocentric distance r by 
r = \JP + Tq — 2/r0 cos ip with Tq = S.bkpc the distance of the Sun to the galactic center 
and ip the direction of the source from the galactic center. AQ represents the solid angle 
at the direction of for a given angular resolution of the the detector. We notice that the 
integration in Eq. Q depends only on the distribution of the dark matter p(r), taken as a 
spherically-averaged form, which is determined by numerical simulation or by observations 
and has no relation to the particle nature of the dark matter. We define this factor as 
'cosmological factor' and the other part in Eq. Q the 'particle factor' which is exclusively 
determined by its particle nature, such as the mass, strength of interaction and so on. The 
factorization of the expression for the annihilation fluxes into a cosmological part and a 
particle part greatly simplifies our discussion. We will discuss the two factors in the next 
sections one by one. 



III. GAMMA RAY FLUX FROM THE GALACTIC CENTER 

In this section we will study the uncertainties from the 'particle factor' of Eq. ((21) by 
calculating the 7-ray flux from dark matter annihilations at the GC. This factor is exclusively 
determined by particle physics and irrelevant to the source of the 7-rays. Several relevant 
concepts will be introduced in this calculation. Since there are observations at the GC we 
will constrain the SUSY parameter space from these observations. 

The Galactic Center is the most extensively studied region for the dark matter anni- 



hilations 



4l| and taken as the most promising source to detect the annihilation products 



since the density cusp at the GC can greatly enhance the annihilation fluxes. However, 



5 



due to the complexity of the GC with the supermassive black hole and many baryonic 
processes there is no convincing signal for the dark matter annihilation even though the 



excess of hig. 
experiments 



n eners 



33, 



y 7-rays beyond the expected background have been detected by several 



3J]. On the contrary, if 7-rays are detected from the subhalos which are 



otherwise completely dark it is a clear signal of the dark matter annihilation. In the follow- 
ing sections we will see that there are other advantages to detect dark matter annihilations 
from the MW substructures. 

We first determine the 'cosmological factor' of the GC in the subsection A and then 
calculate the 'particle factor' in the subsection B by scanning the low energy SUSY parameter 
space and considering the constraints. Finally we give our results in the subsection C. 



A. Cosmological factor 

1. density profile of the MW 

The DM density profile is extensively studied by numerical simulations. However, there 
are still a lot of debates on this subjects, which focus on the slope of the central cusp of 
the profile. It is first given by Navarro, Frenk, and White [42] and supported by recent 
studies 0] that the DM profile of isolated and relaxed halos can be describe by a universal 
form 

pDM{r) = 1—, — rr^ -, — , (3) 

{r/rs){l + r/rsY 

where pg and are the scale density and scale radius respectively. The two free parameters 
of the profile can be determined by the measurements of the virial mass of the halo and 
the concentration parameter determined by simulations. The concentration parameter is 
defined as 

c = ^ , (4) 

where r^i^ is the virial radius of the halo and r_2 is the radius at which the effective loga- 
rithmic slope of the profile is —2, i.e., j;:{r^p{r)) r=r_2 = 0. For the NFW profile we have 
Ts = r_2- The concentration parameter reflects how the DM is concentrated at the center. 
For a larger concentration parameter the DM is more centrally concentrated. 

The NFW profile in Eq. (jHl) has a singularity when r towards zero, p(r) r~^, while its 
slope becomes much steeper at large r, p(r) — > r^^ for r ^ r^. 
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However, Moore et al. gave another form of the DM profile |4J| to fit their numerical 
simulation with an increased resolution 

PDM(r) = , . XI .x/' , I XI , (5) 

which has the same behavior at large radius as the NFW profile while it has a steeper central 
cusps pir) — » r^^'^ for small r than the NFW profile. The index of the central cusp at about 
1.5 is also favored by following higher resolution simulations |. At the present time it 
seems that the different central cuspy behavior are not due to finite numerical resolution of 
the simulations and can not be solved by improving the numerical resolution further. For 
the Moore profile we have = r_2/0.63. 

Anyway, although there is on going debate over which profile is most accurate it is 
reasonable to believe that the NFW and the Moore profiles represent two limiting cases 
between which a realistic description of DM distribution will fall. We will calculate the 
gamma ray fluxes from neutralino annihilation by adopting the two profiles. The actual 
fluxes should fall within the two limiting cases. 

2. core radius 

We notice that both the NFW and the Moore profiles have unphysical singularities at the 
GC which may lead to divergent gamma ray fluxes. Therefore the profiles should have a core 
radius, rcore, within which the DM profile should be kept constant due to the balance between 
the very high annihilation rate and the rate to fill the region by infalling DM particles. The 
time scale of the free fall of the DM particles can be approximately given by!4a| 

while the annihilation time scale is 

{av}n^{r core) 

Taking p about 200 times the critical density and (av) ~ lO^^^cm'^s"^ and applying the 
formulas above we then get rcore ~ 10~^ kpc for the Moore profile and rcore ~ 10~^^ kpc for 
the NFW profile. 

At the position of the Sun, which is about 8.5 kpc away from the the Galactic Center, we 
can probe the GC to a radius as small as 0.15 kpc (0.015 kpc) using an instrument with the 
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FIG. 1: The cosmological factor within the sohd angle Afi as function of the core radius for the 
angular resolution of Ail = 10~^, 10~^ respectively. 

angular resolution of 1° (0.1°), corresponding to the solid angle AQ ^ 10^'^(10^^). Therefore 
we expect the annihilation flux within the solid angle AQ will be enhanced greatly if the 
core radius is smaller than the radius that the instrument can resolve, i.e., rcore < ''"res- In 
Fig. ^ we plot the cosmological factor defined as the integration in Eq. Q as a function of 
the core radius. From the figure we can see that the cosmological factor increases quickly 
with decreasing core radius when it is larger than the rres- For smaller rcore < ''"res, the 
cosmological factor rises slowly by further decreasing the core radius. This behavior can be 
understood as below. The flux is proportional to integration of the density square within 
the solid angle 



$ ~ y p (r)r (ir ~ y ij core) dr + j p^{r)r''^dr 



(8) 



which gives that the flux in Moore profile depends on rcore logarithmically while independent 



of Tcore for the NFW case. 



3. angular resolution 

From Eq. (0) we can see that the cosmological factor also depends on the angular 
resolution of the experimental instruments. We expect a larger cosmological factor for a 
larger angular resolution which probes greater volume of dark matter annihilation (however, 
the significance of detection is decreased due to more background is included, see our later 
discussion). In Fig. we plot the cosmological factor for the solid angular resolution 
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AQ = 10~^ and AQ = 10~^ respectively. For the case when rcore > ''"res the cosmological 
factor for AQ = 10~^ is indeed two orders magnitude larger than that for AQ = 10~^. While 
for smaller rcore the difference between the two cases is smaller. Especially for the Moore 
profile the difference is almost negligible. 

B. Particle factor 

We now turn to the particle factor in Eq. Q. We will mainly work in the frame of MSSM, 
which is a low energy effective description of the fundamental theory at the electroweak scale. 
For comparison, the minimal supergravity (mSUGRA) model is also explored, where the 
soft SUSY breaking parameters can be universally defined at the scale of grand unification 
(GUT). For the R-parity conservative MSSM, the lightest supersymmetric particle (LSP), 
generally the lightest neutralino, is stable and is an ideal candidate of dark matter. 

Even the R-parity conservative MSSM is described by more than one hundred soft su- 
persymmetry breaking parameters. However, for the processes related with dark matter 
production and annihilation, only several parameters are relevant under some simplifying 
assumptions, i.e., the higgsino mass parameter fi, the wino mass parameter M2, the mass 
of the CP-odd Higgs boson m^, the ratio of the Higgs Vacuum expectation values tan (3, 
the scalar quark mass parameter rriq, the scalar lepton mass parameter mj, the trilinear 
soft breaking parameter At and Ab. To determine the low energy spectrum of the SUSY 
particles and coupling vertices, the following assumptions have been made: all the sleptons 
and the squarks have common soft-breaking mass parameters mj and rriq respectively; all 
trilinear parameters are zero except those of the third family; the bino and wino have the 
mass relation. Mi = 5/3 tan^ 6wM2, coming from the unification of the gaugino mass at the 
grand unification scale. 

We perform a numerical random scanning of the 8-dimensional supersymmetric parame- 
ter space using the package DarkSUSY The ranges of the parameters are as following: 
bOGeV < M2, Ma, rriq, mj < 5TeV, 1.1 < tan/3 < 55, -3m^ < At, Ab < Sm^, 
sign(/i) = ±1. The parameter space is constrained by the theoretical consistency require- 
ment, such as the correct vacuum breaking pattern, the neutralino being the LSP and so on. 
The accelerator data constrains the parameter further from the spectrum requirement, the 
invisible Z-boson width, the branching ratio of 6 ^ 57 and so on adopting the 2002 limits 
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of the Particle Data Group j48|. 

Another important constraint comes from cosmology. Combining the recent observation 
data on cosmic microwave background, large scale structure, supernova and data from HST 
Key Project the cosmological parameters are determined quite precisely. Especially, the 
abundance of the cold dark matter is given by 2] ^cdm^^ = 0.113to;Q09. We constrain the 
SUSY parameter space by requiring the rehc abundance of neutralino < Q^h"^ < 0.137, 
where the upper limit corresponds to the 3a upper bound from the cosmological observations. 
When the relic abundance of neutralino is smaller than a minimal value the neutralino 
represents a subdominant dark matter component. We then rescale the galaxy dark matter 
density as p(r) — > ^p(r) with ^ = fli./i^/(^^;^/i^)mm- We take {^x^'^)mm = 0.086, the 3a 
lower bound of the CDM abundance j^. The effect of coannihilation between the fermions 
is taken into account when calculating the relic density numerically. 

Besides exploring the SUSY parameter space at low energy directly under some simplify- 
ing assumptions, there is another popular approach of handling the phenomenology of MSSM 
by assuming a simple SUSY breaking pattern at the GUT scale. One of the simplest models 
in this kind is the minimal supergravity model which has only five free SUSY parameters, 
the gaugino masses mi/2, the sfermion masses rriQ, the trilinear parameter Aq, which are all 
defined universally at the GUT scale, as well as the ratio of the vacuum expectation values 
tan (3 and the sign of the higgino mass parameter fi. The high energy universal relationship 
leads to constraints on the low energy spectrum |49]. To calculate the neutralino annihi- 
lation, we adopt the package ISASUGRA (version 7.69) to calculate the low energy SUSY 
spectrum by numerical solving the renormalization group equations from the GUT scale 
downwards to the weak scale [50|. We also randomly scan the free parameter space consid- 
ering the theoretical and experimental constraints. The ranges of the parameters are given 
as 50GeV < mo,mi/2 < 5TeV, — 3mo < Aq < 3mo, 1.1 < tan/? < 55 and sign(/i) = ±1. 

The 7-rays from the neutralino annihilation arise mainly in the decay of the neutral pions 
produced in the fragmentation processes initiated by tree level final states. The fragmenta- 
tion and decay processes are simulated with Pythia package |5l| incorporated in DarkSUSY. 
We focus our calculation on the continuum 7-rays from the pion decays. 
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1. constraints from cosmic ray observations 



Before presenting the numerical results we discuss the constraints from the cosmic ray 
observations first. High energy and very high en ergy 7-ray emission from the GC have been 



m and HESS 



3J] experiments. However, 



detected by the EGRET |3^, CANGAROO-II 
the present situation is still confused since all these results can not be attributed to a unique 
7-ray source or due to a single emission mechanism. Therefore we take all these results as 
a constraint on the dark matter annihilation, i.e., the predicted 7-rays flux due to the DM 
annihilation should not exceed these observed fluxes. 

The 7-ray spectrum from the EGRET observation 32] at the GC with the angular resolu- 
tion of ~ 1° is well described by a broken power law with a break energy at 1.9GeV . Above 



this energy the photon spectrum is F{E) = 1.6 x 10 ^{E/GeV) 



-3.1 



cm 



^'s ^GeV ^. Assum- 



ing that this spectrum can extend to a quite high energy we get an approximate integrated 
7-ray flux above IGeV using this power spectrum F{E > IGeV) = 7.7 x 10"'' cm~^s~^. 

Both CANGAROO-II and HESS observed 7-rays at the GC from about 200GeV to 
about 10Te\^. The HESS data ^34\ is fitted by a power law spectrum F{E) = 2.5 x 



lO-'^E/TeVy-^^cm-h-'TeV-^ 



cm -l ev ^ The CANGAROO-II |33|, [52 1 gives a quite different 
power law spectrum with a spectral index —4.6. The apparent discrepancy seems indi- 
cate a significant change of the source at lower energies over about one year. However, none 
of the individual experiment observes significant variability of the source. Another possible 



t 



he CANGAROO-II there may 
. The flux detected by HESS 



explanation is that due to larger positional uncertainty o 
be more than one source exist in the direction of the GC I 
is also much lower than that detected by EGRET if extending the spectrum to lower en- 
ergies. Therefore we take a relaxed constraint from the HESS experiment, which has the 
best angular resolution of 0.1°: assuming that the source at the GC can extend to the range 
within an angular resolution of 1°, despite a point-like source can fit the single HESS ob- 
servation welipj]. We then get the integrated flux above lOOGeV using the given spectrum 
F{E > lOOGeV") = 3.4 x 10-^cm~'^s-\ The integrated flux from the CANGAROO-II data 
by extending its spectrum to lower energy, F{E > lOOGeV^) = 2.8 x lO^^cm^^s^^. Anyway, 
at the moment the situation is not clear we take these values as a conservative constraints so 
that we can probe more SUSY parameter space. For further more severe constraints more 
SUSY parameter space will be constrained and our later results can be simply rescaled due 
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FIG. 2: The integrated 7-ray fluxes from neutralino anniliilation at tlie GC above the threshold 
energy of 1, 10, 50, 100 and 1000 GeV respectively for the Moore (upper panel) and the NFW 
(lower panel) profiles as function of the neutralino mass. The fluxes are given within the angular 
resolution of AQ = 10^'^. Each dot in the figure represents a set of low energy SUSY parameters 
which survive all the current limits. The bigger points superposed on each set of dots (corresponding 
to a threshold energy) are given by the mSUGRA model. 



toEq. 0. 
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C. Results 



Combining all these constraints together, from the theoretical consideration to the accel- 
erator experiments and cosmic ray observations, we can give the predicted 7-ray fluxes from 
neutralino annihilations at the GC now. In Fig. |2l we plot the integrated 7-ray fluxes within 
the solid angle AQ = 10~^ as a function of the neutralino mass taking the threshold energy 
as IGeV, lOGeV, 50GeV, lOOGeV and ITeV for both the Moore and the NFW profiles. 
Each dot in the figure corresponds to a model with a set of definite SUSY parameters in 
the 8-dimensional parameter space allowed by all the constraints. The bigger points super- 
posed on each set of dots are given by the mSUGRA model. We notice that the mSUGRA 
model tends to give large annihilation fluxes, while small neutralino mass which is hard to 
be greater than SOOGeV. The scatter of the points represents the uncertainty coming from 
the unknown soft SUSY breaking parameters. 

It is worth giving some comments here. First, the core radius is taken as lO"*^ kpc for 
both the Moore and the NFW profiles. From Fig. Qwe can see that the flux in the Moore 
profile decrease only by a factor of 2 when taking rcore from 10~^ kpc to 10~^ kpc. For the 
case of NFW profile the predicted flux has no change even taking r^ore as large as 0.1 kpc. 
Second, in the case of Moore profile the SUSY parameter space has been constrained by the 
cosmic ray observations at the GC. If more stringent constraints are adopted all the later 
results according to the Moore profile will be rescaled. However, there is no constraints 
for the case of the NFW profile. Third, except the region near the threshold energy the 
uncertainties of the predicted 7-ray fluxes are within two orders of magnitude even scanning 
the 8-dimensional parameter space which spans quite a large range, as given in the last 
subsection. The most strong constraint comes from the requirement by cosmology, i.e., 
requiring < ^cdm^^ < 0.137. For s-wave annihilation the rate (erf )ann is closely related 
with the initial value of (af )ini at the decoupling epoch which determines the relic density 
of neutralino. Forth, the neutralino mass above 1 TeV is difficult to achieve after taking all 
the constraints into account. We can realize a model with the neutralino mass at most as 
heavy as about 5 TeV to satisfy all the constraints after we relax the range of all low energy 
soft mass parameters as large as 50TeV. Therefore, a neutralino as heavy as 18 TeV to 



explain the spectrum observed by HESS experiment 



53, 



54 , n5| as neutralino annihilation 



will be very hard to achieve in the framework of MSSM. 
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For discussions of the substructure emission in the next sections we will fix the 'parti- 
cle factor' by taking an optimistic set of SUSY parameters. The 'particle factor' is such 
taken that the integrated 7-ray fiuxes above lOOGeV^ from the GC for the Moore profile is 
F{E > lOOGeV) = lO^^cm^^s^^ , which is near the maximal value of the SUSY prediction. 
This is just a normalization of the 7-ray fiux and only the relative magnitude of the fiuxes 
between the GC and the substructure is relevant, since the overall magnitude can be rescaled 
according to different particle factors. 

We can summarize this section here. By scanning the SUSY parameter space and after 
taking all the constraints into account we give the scatter of the integrated 7-ray fiuxes from 
neutralino annihilations at the GC in Fig. |21 for different threshold energies and for both 
the Moore and the NFW profiles. From the factorization of the expression in Eq. we 
know this figure is suitable to any other source by multiplying each point by a global factor 
which represents the difference of the 'cosmological factor' from that at the GC. Therefore 
the figure shows the general uncertainties of the particle factor. The uncertainty is quite 
small considering the huge volume of the free parameter space. This is due to the fact that 
the annihilation process is closely related with the process of dark matter freezing out which 
determines its relic density. 

IV. GAMMA RAYS FROM THE SUBHALOS 
A. realization of MW with substructures 

To predict the 7-ray fiux by neutralino annihilation from the subhalos we need to know 
the distribution of subhalos in the MW and the density profile within each subhalo. 

The properties of subhalos are determined via the competition between accretion and 
destruction due to tidal force and dynamical friction. N-body simulation and semi-analytical 
methods have been extensively used to investigate the spatial distribution and mass function 
of substructures in the host halo. According to the extensive studies it is now believed 
that the radial distribution of substructures is generally shallower than density profile of 
the smooth background. The reason for the anti-bias of the substructure number density 
relative to the smooth distribution is due to the tidal disruption of substructures which is 
most effective near the galactic center. This conclusion does not depend on the numerical 
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resolution as confirmed in Ref. [56[ by adopting a wide range of mass and force resolutions. It 
is shown that the relative number density of subhalos can be approximated by an isothermal 
profile with a core 

n(r) = 272^(1 + {r/rHfr^ , (9) 

where nu is the relative number density at the scale radius th- The average core radius for 
the distribution of galaxy subhalos is about 0.14 times the halo virial radius, th = 0.1 4?",,; 
The core radius is a smaller fraction of the virial radius than that of cluster subhalos jSy], 
since galaxy forms earlier and is more centrally concentrated. The spatial distribution given 
above agrees well with that in another recent simulation by Gao et al. js^]. 

At large radius n{r) in Eq. © goes as which represents the anti-bias since the DM 
profile goes as at large radius for both the NFW and the Moore profiles. It is worth 
mentioning that most previous works calculating 7-rays from substructures adopting a fitted 



formula given in Ref. 58|, n(r) oc (1 + (r/r/^)^)"^'^, which follows the DM profile and does 
not reflect the anti-bias of the substructure distribution. 

Simulations show that the differential mass function of substructures has an approximate 
power law distribution, dn/dm ~ m~°, with no dependence or slight dependence on the host 
halo mass. Most studies point out that the cluster and galaxy have similar substructure mass 
function although the clusters form much later than galaxies in the hierarchical structure 
formation scenario. It seems that the tidal effects change the mass distribution function self- 
si,„,.a.. In Ref. H bo*K *e dust, and galaxy sub.t.nct.e cumulative ,„as. fuuetious 
are found to be an m^^ power law, risub(^sub > m) (x m~^, with no dependence on the mass 
of the parent halo. A slight difference is found in a .cent .„.,at.n.; Gao.. .RtHat 
the cluster substructure is more abundant than galaxy substructure since the cluster forms 
later and more substructures have survived the tidal disruption. The mass function for both 
scales are well fitted by dn/dm oc m^^'^. A non-universal form of the mass function is found 

n 

in Ref. |6']: the power index of the cumulative mass function is —1 for mgub > lO^^/i ^Mq 
while it changes to —0.7 for msub < IO^^H^^Mq. There is an advantage of a power law form 
for the differential mass function shallower than m~^: the fraction of the total mass enclosed 

raction 



in subhalos is then insensitive to the mass of the minimal subhalo we take. The mass : 
of subhalos estimated in the literature is around between 5 percent to 20 percent 
In this work we will always take the mass fraction of substructures as 10 percent. 

Putting all these arguments together, we get the probability of a substructure with mass 
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m at the position r to the galactic center 



n(m, r) = rio j (1 + {r/rnf)-' , (10) 

where Mvir is the virial mass of the MW, no is the normahzation factor determined by 
requiring the total mass of substructures converges to 10 percent of the MW virial mass, 
Mvir- A population of substructures within the virial radius of the MW are then realized 
statistically due to Eq. ()10|). The mass of the substructures are taken randomly between 
Mjiin = IO^Mq, which is the lowest substructure mass the present simulations can resolve 
[Sd|, and the maximal mass Mmax- The maximal mass of substructures is taken to be 
0.01 Mvir since the MW halo does not show recent mergers of satellites with masses larger 
then ~ 2 X 10^^ Mq. It will be shown that the 7-ray flux is quite insensitive to the minimum 



subhalo mass since the flux from a single subhalo scales as its mass j30l l61|. 

We notice that the number density of substructures is largest at the GC due to Eq. (|Tn]l 
which, however, is in conflict with the fact that most substructures near the GC are destroyed 
by the strong tidal effect. The underestimate of the tidal effect near the GC in Eq. (jlOj) is 
due to the finite resolution of the N-body simulations and the formula is an extrapolation of 
the subhalo distribution to smaller radius. The global tides from the host halos strip outer 
parts of the substructures and result in the substructure disruption or at least significant 
amount of substructure mass loss. We take the tidal effects into account under the "tidal 
approximation" , which assumes that all mass beyond a suitably defined tidal radius is lost 
in a single orbit while keep its density profile inside the tidal radius intact. 

The tidal radius is defined as the radius of the substructure at which the tidal forces of 
the host exceeds the self gravity of the substructure. Assuming that both the host and the 
substructure gravitational potential are given by point masses and considering the centr ifug al 
force experienced by the substructure the tidal radius at the Jacobi limit is given by j6^ 

'''' = (^) ^ ' ^''^ 
where is the distance of the substructure to the GC. The substructures with rtid < will 
be disrupted completely and be discarded in our realization of the substructure population. 
The mass of a substructure is also recalculated by subtracting the mass beyond the tidal 
radius. Therefore the final radial distribution of substructures near the GC is somewhat 
different from that given in Eq. which is shown in Fig. El Indeed the substructures 
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FIG. 3: Number of substructures along the radius of the MW for the Moore and the NFW profiles. 

near the GC are disrupted completely after we take the tidal effects into account. The 
substructures with NFW profile can exist more near the GC than the Moore profile. This 
is because that the NFW profile is more centrally concentrated with smaller r^. 

B. Concentration parameter 

Once the population of the substructures are determined we need to determine the dark 
matter density profile in each substructure in order to calculate the annihilation of neu- 
tralinos. The density profile in the substructure, either the NFW or the Moore profile, is 
characterized by two scale parameters ps and r^. The two parameters can be determined us- 
ing two different methods. In the first method they are determined by one virial parameter, 
the virial mass (or equivalently the virial radius, or virial velocity), and the concentration 
parameter, which relates the virial and the inner scale radius. 

N-body simulation shows that the concentration of the substructure is strongly correlated 
with the formation epoch of the substructure. At an epoch of redshift Zc a typical collapsing 
mass M^{zc) is defined by (t[M*(2;)] = 5sc{z)i where the o"[M*(2)] is the linear rms density 
fluctuation on the comoving scale encompassing a mass M,,, 5sc is the critical overdensity 
for collapsing at the spherical collapse model. The collapsing mass M^, is determined by 
the primordial linear power spectrum of the fluctuations and the known cosmology, which 
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FIG. 4: Concentration parameter as a function of the virial mass of the halo calculated according 
to the Bullock modelfSj]. The model parameters are taken as -F = 0.015 and K = 4.4. The 
cosmology parameters are taken as = 0.3, Oa = 0.7, Qsh^ = 0.02, h = 0.7, ag = 0.9 with 
three generations of massless neutrinos. 

determines the evolution of the fluctuations. In a semi- analytic model Bullock et al. 0] relate 
the typical collapsing mass to a fixed fraction of the virial mass of a halo M^:{zc) = FM^i,.. 
The concentration parameter of a halo with virial mass Mvir at redshift z is then determined 
as Cvir(Mvii., 2;) = K^^^. Both F and K are constants to fit the numerical simulations. 
Once F and K are determined the concentration of a halo is completely determined by the 
cosmology in hand. We notice that a smaller Mvir corresponds to a smaller collapsing mass 
and early collapsing epoch when the Universe is denser and therefore a larger concentration 
parameter. Therefore smaller subhalos are more closely concentrated. 

In Fig. 13] we plot the concentration parameter at 2; = as a function of the virial mass 
of a halo according to the Bullock model 6^ . In the calculation, we have taken a standard 
scale invariant primordial spectrum of the fluctuation with the cosmology parameters taken 
as Qm = 0.3, Qa = 0.7, fi^/i^ = 0.02, h = 0.7, as = 0.9 and three generations of massless 
neutrinos. The model parameters are taken as F = 0.015 and K = 4.4. The scatter of 
the concentration parameters for a given halo mass is log-normal with la deviation around 
the mean as A(logioCvir) = 0.14. In Fig. |3] both the median and the ±la values of the 
concentration parameters are plotted. 
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The Bullock model reproduces the concentration parameters quite well from the N-body 
simulations. From Fig. |3] the experiential formula is confirmed again that c^h oc M^f . 
We expect that this exponential relation of the concentration parameter and virial mass 
for subhalos should be very well followed, since subhalo forms early at the epoch when the 
Universe is dominated by matter with approximate power-law power spectrum of fluctua- 
tions. Besides the Bullock model we have also adopted other recent simulation results in 
the literature. We use the experiential relation between the concentration parameters and 
virial mass and fit the parameter f3 from these simulation results. The density profile of the 
substructure and furthermore the 7-ray flux from the substructure are then calculated. We 
find the concentration parameter is the most sensitive parameter in determining the 7-ray 
flux. Different behavior of the concentration leads to different predictions of the 7-ray fluxes. 

The second method to determine the profile parameters, as given in requires each 
substructure produced following Eq. (fTUj) is compact enough to resist the tidal stripping. 
Since the distribution of substructures in Eq. (|T0|) is given by simulations at z = the 
second method reflect the simulation results faithfully in the statistical meaning. Since the 
simulation can not resolve the region near the GC we have given a cutoff at r = 10 kpc from 
the GC requiring no subhalo exits inside the cutoff radius. The conditions to determine the 
parameters are given by|6^]: /2vir = -Rtid, I^"" p{r)dV = Msub and phost{rd) = Pd(-Rvir)- The 
last condition requires that the density of the subhalos at R^r should equal the local density 
of the host halo at the position of the subhalo ra- 

C. Results 

We can now present the results of the 7-ray fluxes from the neutralino annihilation in 
the MW substructures. The fluxes are averaged within the angular resolution of 1° in the 
following way: if a subhalo is too small or too far from the detector that its angular size is 
smaller than the angular resolution we calculate the annihilation flux of the whole subhalo; 
on the contrary for these subhalos the detector can resolve we calculated the flux within 
the 1° angular resolution assuming that we are aiming at the center of the subhalo. The 
annihilation core radius in the following calculations is taken as Vcut = 10~^kpc, which is 
quite a conservative value and affects the the results little according to Fig. [T] 

In Fig. Elwe plot the integrated 7-ray fluxes above the threshold energy of 100 GeV as a 
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FIG. 5: The integrated 7-ray fluxes above lOOGeV as function of the angles of the sources relative 
to the GC. In the upper panel we adopt the concentration parameters of the subhalos given in [g^ 
while in the lower panel we adopt that in ^5}. Both contributions from the smooth component of 
DM and from subhalos are plotted. 



function of tp, the angle of the source to the direction of the GC. In the figure no information 
is shown about the other direction around the axis along the Sun and the GC. Only the 
maximal 7-ray flux at this direction is chosen for each ip and plotted in the flgure. Both 
the 7-rays annihilated from the smooth dark halo and that by adding the emission from 
the subhalos and the smooth component are shown for the Moore and the NFW proflles. 
For the smooth component the 7-ray flux decreases rapidly when departing from the GC 
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while for the subhalos the contribution seems almost isotropic to different directions. This 
is due to the fact that the Sun is only 8.5 kpc away from the GC. In the upper panel 
we calculate the flux adopting the Bullock modeljcsl and in the lower panel we adopt the 
concentration parameter according to the simulation by Reed et al. pa]. Generally the 
Moore profile predicts greater 7-ray fluxes. At the GC the 7-ray flux for the Moore profile 
is about 300 times higher than that for the NFW profile. However, for the subhalos the 
difference between the two profiles is only about one order of magnitude. Therefore the 
uncertainty of the predicted 7-ray fluxes from the subhalos due to the cosmological factor 
is much smaller than that from the GC. The figure can be used to other threshold energies 
by a global shift of the particle factor from the Fig. |21 

It should be noted that Fig. El illustrates the 7-ray fluxes from the subhalos due to a 
statistical realization of our Galaxy. None of the positions of peaks are predicted exactly and 
the maximal flux may have a large fluctuation: it is possible that a subhalo is accidentally 
located near the solar system. Therefore we try to calculate the statistically averaged fluxes 
by realizing one hundred MW sized halos and count the average number of subhalos emitting 
7-rays above any intensities. 

Fig. ini gives the cumulative number of subhalos emitting 7-rays with intensity above the 
integrated flux It should be noted that the results, also for other figures hereafter, are 
given within the zenith angle of 60°, which is the maximal angle a ground array can possibly 
probe, instead of the whole sky. In the upper panel we plot the results for the Moore profile 
while the lower panel is for the NFW profile. The curves are given by calculating the density 
profile of subhalos according to different author's simulation results, where 'cl' denotes the 
simulation of Ref. Q; 'c3' of Ref. Q; 'c4' of Ref. H for the ACDM model with as = 0.9; 

c6' uses the median Cmr — M^r relation for distinct halos of the Bullock model given in Ref. 

6^, while 'c7' and 'c8' take the upper 2(T and la values of the same model respectively; 
'c2' adopts the simulation results given in dense matter environment of the same reference 
and extends the relation to small subhalos; since this relation gives very large concentration 
parameters for small subhalos we have cut the concentration parameter arbitrarily if c > 100 
in 'c5'. From this figure we can easily read the number of the expected detectable subhalos 
if the sensitivity of a detector is given (with same threshold energy and angular resolution 
adopted here). 

We find the cumulative number of subhalos for large fluxes can be well fitted by an inverse 
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FIG. 6: The cumulative number of subhalos as function of the integrated 7-ray fluxes n(> 
for the Moore profile (upper panel) and the NFW profile (lower panel). The results are given 
within the zenith angle of 60°. The curves represent the results according to different simulations 
as explained in the text. These curves give the number of subhalos which emit 7-rays with the 
integrated flux above <I>^. 



power law, n{> $^) = uq/^^, as shown in Fig. [7| The deviation of the calculated curves 
from the fitted lines at the end of largest fluxes is due to the fact that we do not have 
enough statistics here, while the deviation at the end of lowest fluxes may be due to the cut 
of the minimal mass of subhalos. In the table H] we give the fitted constant uq for different 
curves. The Moore profile predicts about 3 times more detectable subhalos than the NFW 
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FIG. 7: The cumulative number of subhalos as function of the integrated 7-ray fluxes n(> 
for the Moore profile (upper panel) and the NFW profile (lower panel). The results are given 
within the zenith angle of 60°. The lines are the corresponding fits of the curves for high 7-ray 
fluxes. From these fits we can easily read out the number of subhalos which emit 7-rays with the 
integrated flux above <I>^. 



profile averagely. Therefore the 7-ray fiuxes from the subhalos are not so sensitive to the 
dark matter profiles as these from the GC, as shown in the Fig. El In Fig. IHlwe show how 
the cumulative number of subhalos changes as the value we take for the minimal mass of 
subhalo. We illustrate the result in the model 'cl'. The dependence on the minimal mass 
of subhalo is very weak, which can be understood from some simple scaling arguments js^. 
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FIG. 8: The cumulative number of subhalos as function of the integrated 7-ray fluxes n(> for 
taking different cuts of the minimal mass of the subhalos. The results are given within the zenith 
angle of 60°. The lines are the corresponding fits of the curves for high 7-ray fluxes. The results 
are quite insensitive to the cut of the minimal subhalo mass. 





cl 


c2 


c3 


c4 


c5 


c6 


c7 


c8 


Moore 
NFW 


1.5 X 10"^^ 
4.5 X 10-12 


2 X 10"^ 
7 X 10-10 


6 X 10^11 
2 X 10-11 


8 X 10-12 
3.5 X 10-12 


3 X 10-1° 
9 X 10-11 


2 X 10-11 
6 X 10-12 


8 X 10-11 
3 X 10-11 


4 X 10-11 
1.5 X 10-11 



TABLE I: Values of the constant no in the fit n(> <I>^) = tiq/^^i to the curves in Fig. El for the 
Moore and the NFW profiles. 



It is worth noting that taking smaller minimal subhalo mass indeed makes the curves closer 
to the fitted lines at the end of smallest fluxes. 

Finally, in Fig. IHlwe give the results according to the second method of determining the 
concentration parameter. Both the NFW and Moore profiles predicts much larger fluxes 
compared with that from the model 'c6' using the Moore profile, while smaller than that of 
the 'c2' model. 

In summary, we calculated the 7-ray fluxes from the MW substructures in this section. 
There are extensive studies on the evolution and distribution of substructures by numerical 
simulations. We integrate the recent simulation results in our calculation. The different sim- 
ulation results give uncertainties in predicting the 7-ray fluxes. The concentration parameter 
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FIG. 9: The cumulative number of subhalos as function of the integrated 7-ray fluxes n(> $^). 
The results are given within the zenith angle of 60°. The curves are calculated according to the 
second method to determine the profile parameters as described in the text. The corresponding 
results from the models 'c2' and 'c6' are also shown. 

is the most sensitive parameter in determining the 7-ray fluxes. We give the statistical aver- 
age number of subhalos with 7-ray intensity above some values. The result is well fitted by 
an inverse power law. From these results the number of detectable subhalos is easily read 
out if the sensitivity of a detector is known. We will discuss the detectability of the 7-ray 
signals from the MW subhalos in the next section. 

V. DETECTABILITY 

The detectability of a signal is defined by the ratio of the signal events to the fluctuation of 
the background. Since the background follows the Poisson statistics, its fluctuation has the 
amplitude proportional to \/Nb- The significance of the detection is quantified by a = 



where eAf^ = 0.68 is the fraction of signal events within the angular resolution of the in- 
strument and the integration is for the energies above the threshold energy Eth-, within the 
angular resolution of the instrument AVL and for the observational time. The effective area 



The signal events are given by 




(12) 
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Aeff is a function of energy and zenith angle. The (f){E) is the flux of 7-rays from DM 
annihilation. 

The corresponding expression for the background is similar to Eq. (|12|) . The background 
includes contributions from the hadronic and electronic comic-rays and the Galactic and 
extragalactic 7-ray emission. We have adopted the expressions as 

(f)h{E) = lA9E-^'^*cm-^s-^sr-^GeV^^ (13) 

for the hadronic contribution Gsl ]. 

(j)eiE) = 6.9 X IQ-"^ E-^-^cm-'^s-\sr-^GeV-^ (14) 

for the electronic contribution 6^, 

(/.extra-7(^) = 1-38 X 10'^ E'^'^ 0711'^ s^Ksr-^ GeV'^ (15) 

for the extragalactic 7-ray emission extrapolated from EGRET data at low energies 13] and 

0gaiac-7(^) = No{l, b) X lO^^E^^'^cm^^^-^sr'^G'el/"^ (16) 



for the Galactic 7-ray emission, also extrapolated from the EGRET data at low energies|71l|. 
with Nq the normalization factor depending on galactic coordinates (/, b). The A^o is modeled 
using EGRET data at 1 GeV [zil 

iVo(^ b) = , ^^'^ + 0.5 (17) 

^1 + (//35)' VI + (6/(1.1 + \l\ 0.022))' 

for |/| > 30° and 

iVo(/, b) = + 0.5 (18) 

v/l + (//35)' Vl + (6/1.8)' 

for |Z| < 30°, with the longitude / and the latitude 6 varying in the intervals [—180°, 180°] 
and [-90°, 90°], respectively. 

Since most background comes from the hadronic cosmic rays, the hadron-photon identi- 
fication efficiency is an important factor to reduce the physical background. For a satellite 
borne experiments, such as GLAST jssj], an identification effici ency of charged particles 
as high as 99.997% can be assumed, while 90% for the photons|72]. However, the effec- 
tive area of this kind of experiments is limited by the size of the satellite and has the 
order of Aeff ~ Im^. The atmospheric Cerenkov telescopes (ACT), such as VERITAS js^, 
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FIG. 10: The 3a sensitivity curves of different detectors as functions of the threshold energy of 
the 7-rays. The curves are calculated according to the parameters given in the text. 
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381 1 , can have very large effective area with an identification efficiency 
of 99% for both the hadronic and the electromagnetic primary particles. However, the ACTs 
have a small field of view (~ 5°) with a duty cycle of about 10%. Therefore the ACTs are 
suitable for the observation of point sources and can not do the blind search. 

The ground-based extensive air shower (EAS) arrays, such as ARGO^], MILAGRoH 
and the next generation all-sky high energy 7-ray telescope HAWcfl, 

have complementary 

properties to the satellite borne experiments and the ACTs. They also have large effective 
areas and at the same time they have large filed of view (~ 60°) and a duty cycle of about 
100%. However, the EAS arrays have low hadron-photon identification efficiency. For the 
ARGO we assume no discrimination between the hadron and the photon, while for the 
HAWC the hadron-photon discrimination can improve the significance of the detection by 
a quality factor of 1.6 |40|. 

In Fig. ^1 we show the 3a sensitivity curves of different detectors as function of the 
threshold energy. We have assumed the angular resolution as 0.1° for GLAST and VERITAS 
and the angular resolution as 1° for ARGO and HAWC. For GLAST, ARGO and HAWC 
we use 1 year of data taking, while for VERITAS 250 hours observation pointing to a 
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source is assumed. We make rough approximation of energy independent effective area, 
lO^m^, 3 X lO^m^ and 4 x lO^m^ for ARGO, HAWC and VERITAS respectively. The diffuse 
extragalactic 7-ray background from Eq. (jl5|) is also plotted in the figure for comparison. 

The detectability of GLAST and VERITAS for the 7-rays from the Galactic subhalos 
has been carefully studied in the Ref. 30, It is shown that the significance of the 
detection of the annihilated 7-rays at GLAST is maximized for ^ AOGeV, since GLAST 
has low energy threshold and thus introduces large background and the 7-ray flux decreases 



when increasing the neutralino mass. For neutralino heavier than lOOGeV^ GLAST has little 
chance to detect the signal js^. We can get similar conclusion from Fig. |21that for small 
threshold energy the maximal flux is at the lower end of the neutralino mass while as the 
threshold energy becomes higher the position of the maximal flux moves to a neutralino mass 
of several hundred GeV. VERITAS, with a higher threshold energy, is sensitive enough to 
detect the 7-rays annihilated from heavy neutralinos. However, its small field of view limits 
its ability to do such observation j^^. 

However, it is possible to do the observation by the EAS array according to our calcula- 
tions. Taking the threshold energy as lOOGeV, ARGO and HAWC require the 7-ray fluxes 
as high as 8 X 10^^^cm~^s~^ and 1 x 10~^^cm~^s~^, respectively, to achieve a 3a significance 
detection. Comparing with the Table I, we get that ARGO can detect 0.1 to as many as 
about 25 subhalos with a Moore profile and from 0.04 to about 9 subhalos with a NEW 
profile for one year data taking. For HAWC, since its sensitivity is about 8 times higher 
than that of ARGO it can also detect more subhalos by 8 times than that of ARGO. For 
example, for the Bullock model and its 2a upper limit, ARGO has the ability to detect 1/4 
and 1 subhalos respectively for the Moore profile and about 0.1 to 0.4 subhalos respectively 
for the NFW profile. While for HAWC even with NEW profile and the median concentration 
parameter of the Bullock model it can detect about 1 subhalo for one year observation. 

Here we can conclude that for the neutralino dark matter as light as about lOOGeV 
GLAST is most suitable to observe the annihilated 7-rays from the Galactic subhalos. For 
the heavier neutralino dark matter ~ 500GeV (< ITeV) ARGO/HAWC can be a viable 
complementary. Especially, if HAWC is built it has great potential to do the observation. 
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VI. DISCUSSION AND CONCLUSIONS 



In this work we calculated the 7-ray fluxes produced by the dark matter annihilation in 
the Galactic subhalos and discussed the detectability of such signals by different types of 
detectors according the nature of the dark matter. 

We explored the low energy parameter space of the MSSM and studied the uncertainties 
from the particle physics in predicting the 7-ray fluxes. Uncertainties from the astrophysics 
are also carefully studied, where we find the most sensitive parameter in predicting the 7-ray 
fluxes is the concentration parameter of subhalos. At the moment there exist discrepancies 
according to different author's numerical simulations. We present the results according to 
several recent simulation results. 

Assuming optimistic SUSY parameters the 7-rays from subhalos may be detected by 
satellite borne experiment, such as GLAST, which has large field of view and small effective 
area if the 7-ray flux is large enough (when the neutralino is light ^ lOOGel^). On 
the contrary, when the neutralino mass is large (~ 500GeV) the 7-ray flux is reduced 
and only ground based experiments with large effective area and large field of view, such 
as ARGO/HAWC, can do the job. We calculated the statistic average numbers of the 
detectable subhalos at these detectors. 

If such an annihilation signal is indeed detected in the future it will not only indicate the 
weakly interaction between the dark matter particles and further implicate the nature of new 
physics beyond the SM but also tell us a lot about the nature of the subhalos: they must have 
a cuspy profile with the Moore or the NFW form or some form between them and the CDM 
scenario is favored without power suppression at the subgalactic scale. However, from this 
single measurement neither the SUSY parameter space nor the subhalos distribution and its 
profile can be actually determined, since both sides still have large uncertainties. Anyway, 
the indirect search of dark matter provides valuable complementary to both the collider 
study of particle physics and the more precise simulations of the dark matter evolution. 

Finally we want to stress again the advantages of search for the dark matter annihilation 
from the MW subhalos. First, subhalos produce clean annihilation signals, as we have 
explained before. The annihilation radiation from the GC is heavily contaminated by the 
baryonic processes. Furthermore, the density profile near the GC is complicated due to the 
existence of baryonic matter. For example, the SMBH can either steepen or flatten the 
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slope of the DM profile at the innermost center of the halo[3|. For subhalos, their profile 
may simply follow the simulation results. Second, the small subhalos form earlier and have 
larger concentration parameter, which leads to relatively greater annihilation fluxes. Third, 
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the DM profile may be not universal, as shown in the simulation given in Ref. 
Smaller subhalos have steeper central cusp. In this case, from Fig. El taking the GC the 
NFW profile and the subhalos the Moore profile, the 7-ray fluxes from the subhalos may 
even be greater than that from the GC. Forth, according to the hierarchical formation of 
structures in the CDM scenario we expect that subhalos should contain their own smaller 
sub-subhalos, which can further enhance the annihilation flux. The sub-subhalos have been 
observed in the numerical simulations, such as in the Ref. 771. Finally, the environmental 
trend seems to make the subhalos more concentrated [63j . However, the effects need further 
studies by more precise simulations. 
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